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Introduction 


With  very  few  exceptions  the  modelling  of  the  compressive  response 
of  fiber-reinforced  composites  is  based  upon  the  proposition  that  fiber  and 
matrix  regions  can  be  considered  as  layered  plates.  In  addition,  the 
deformation  of  the  soft  matrix  material  is  assumed  to  vary  linearly  between 
the  boundary  values  imposed  upon  the  matrix  layer  by  the  neighboring 
fibrous  plates. 

In  this  article  it  is  shown  that  the  latter  assumption  yields  inacceptable 
deformation  fields  within  the  matrix  when  employed  in  conjunction  with 
realistic  fiber  geometries,  such  as  hexagonal  or  square  arrays  of  circular  fibers. 
This  observation  casts  serious  doubts  about  the  validity  of  the  basic  premise 
employed  in  the  modelling  of  the  compressive  respoitse  of  composites. 

Review-of  Current  Modelling 

Let  X,  Y,  z  denote  rectilinear  Cartesian  coordinates  and  consider  a 
composite  material  consisting  of  periodic  layers  of  fiber  and  matrix  materials, 
with  a  cross  section  shown  in  Figure  1. 

Let  N  denote  a  compressive  load  acting  parallel  to  the  z  axis  and 
assume  that  the  composite  deforms  in  a  "shear  mode"*'^  as  shown  in  Figure  2. 
Denote  by  Wrj  and  Wf  the  components  of  displacement  in  the  direction  of  z 
in  the  matrix  and  fiber  layers,  respectively.  The  basic  assumptions  employed 
in  the  modelling  of  the  compressive  response  of  composite  are  expressed  by 

w„(x,z)  =  W,Kz) 

and 

Um(X,z)  =  Uf{c  +  h,z) 

Equation  (2)  derives  from  the  assumption  that  the  transverse 

translation  of  fiber  and  matrix  regions  is  independent  of  X. 

An  essential  simplification  occurs  when  the  fiber  is  assumed  to  rotate 
rigidly  about  its  mid-plane,  in  which  case  we  have 

Wf(-c,z)  =  - Wf{c,z)  (3) 

A  most  significant  consequence  of  expressions  (1),  (2)  and  (3)  is  that  the 
displacements  in  the  matrix  region  take  the  forms 

W„(X,z)=^^i^X  =  F(z)X  (la) 

and 

U„(X,z)  =  P(z)  (2a) 


(1) 

(2) 
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Thereby,  the  shear  strain  in  the  matrix  depends  on  the  coordinate  z 
alone,  namely 


YS  =  F(z)  +  P'(z)  =  Q(z) 


(4) 


Obviously,  we  also  have  YyS  =  0. 

It  will  be  shown  in  the  next  section  that  displacements  of  the  forms 
given  in  equation  (1),  or  (la),  are  inacceptable  for  realistic  fiber  geometries 
because  they  yield  singular  values  for  Those  singular  values  give  rise  to 

unbounded  terms  in  the  strain-energy  of  the  matrix,  and  therefore  render  the 
layered  geometry  unsuitable  for  modelling  the  compressive  response  of  fiber- 
reinforced  composites. 

Hexagonal  Fiber  Arrays 

Consider  the  hexagonal  array  of  circular  fibers  with  a  cross  section 
shown  in  Figure  3.  Let  all  the  rows  of  fibers  parallel  to  the  Y  axis  undergo 
rigid  rotations  which  are  analogous  to  the  rotations  of  the  fiber  layers  shown 
in  Figure  2. 

To  focus  ideas  consider  the  representative  cross-sectional  area  depicted 
in  Figure  4,  and  refer  to  Cartesian  coordines  x,  y  translated  to  the  center  of  that 
area. 

Let  Cl  and  C2  denote  the  arcs  which  separate  the  fibers  and  matrix 
regions  in  the  representative  area,  and  divide  the  matrix  into  three  sub- 

regions  A,  B,  and  A  .  Also,  let  a  denote  the  radius  of  the  fibers  centered  at 

and  (-^b  ,  -  ,  as  shown  in  Figure  4. 


We  have 

(x-^b)^+jy-tj^  =  a2  (5) 

hence 

x{Ci)  =  J^b- Y 

On  C2:  (x  +  +  (y  +  ^)^  =  (6) 

hence 

x(C2)  =  -^b  +  ya2-(y+tj^  (6a) 
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In  analogy  with  the  circumstance  shown  in  Figure  2  we  have 


(7) 


However,  we  now  also  have 


(8) 


In  addition,  in  view  of  the  presumed  rigid  rotations  of  the  fibers  we  obtain 
OnQ: 

w„  =  W,  =  W,  .  F  (z)  (x  (C,)  -  fib)  -  -  F  (2)^  a2-(y-b)^  (9a) 

while  on  C2: 

Wm  «  Wf  =  W2  =  F(2)  (x(C2)  +  ^b)  =  F(z)  Y  a2  - (y  +  (9b) 

Assume  that  the  matrix  in  sub-regions  A,  B,  and  a'  deforms  in 
accordance  with  expressions  (1)  and  (la).  In  this  case  we  obtain  the  following 
expressions  for 


In  sub-region  A: 


Wm  = 


=  Wi  ^^^  =  ~F(z)(x-»-^b)-=^ 
X(Q)  +  J5:b  '  2  /V3b- 


f(y) 
f(y) 


In  sub-region  B: 


(10a) 


=  F(2) 


g(y)-(f(y)  +  g(y)) 


x-g(y)-»-^b 
V3b-(f(y)  +  g(y)). 


(10b) 
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In  sub-region  A*: 


Wm  =  W2 


^h-4C:d  'V3b-g(y) 


(10c) 


In  equations  (9):  f  (y)  =  and  g(y)  =  ya2-(y  +  t]^  . 

From  equations  (10)  it  is  noted  that,  in  clear  contradiction  with 
expressions  (1)  or  (la),  a  linear  variation  of  with  x  in  hexagonally  arrayed 
circular  fibers  must  take  the  form 

W™=F(z)(R(y)x  +  r(y))  (11) 


Instead  of  =  F(z)x. 


The  qualitative  distinction  between  the  displacement  fields  which 
correspond  to  expressions  (1)  and  (10)  is  exhibited  in  Figure  5. 

By  hypothesis  V,„  =  0,  therfore  upon  differentiation  of  expressions  (10) 
with  respect  to  y,  one  obtains  the  following  expressions  for  : 


In  sub-region  A: 


ay 


y -b. 

=  F(z)V3b(x  +  V3b) - 2 - 

f(y)(y3b-f(y))^ 


(12a) 


In  sub-region  B: 


ay 


V3bF(z) 

[lf3b-(f(y)  +  g(y)jl^ 


(f(y)  +  x-fib) 


+ 


b  +  X 


-g(y))} 


(12b) 
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Since  f|^-a|  =  g|a- ^1  =  0,  expression  (12b)  yields  unbounded  values, 
singular  of  power  1/2,  along  the  lines  X*  t(a-^|  which  separate  sub-region 


B  form  subregions  A  and  A.  The  above  singular  terms  result  in  a  non- 
intergrable  (power  -1)  singular  expression  for  the  strain-energy  within  the 
matrix.  Consequently,  the  basic  premises  employed  in  the  "layered 
formulation"  in  equations  (1)  or  (la)  fail  to  provide  a  meaningful 
approximation  to  the  response  of  fiber-reinforced  composites. 


It  should  be  noted  that  Figures  4  and  5  correspond  to  the  case  a  > 


2 


/ 


namely  for  fiber  volume  fractions  Vf  >  0.25. 

The  case  of  Vf  <  0.25  corresponds  to  the  configuration  shown  in  Figure 

6,  in  which  case  the  expressions  for  in  sub-regions  A  and  A*  are  still 
given  by  equations  (10a)  and  (10c)  while  in  region  B  =  0.  In  this  case  we 


again  obtain  singular  values  for 


awn, 

ay 


along 


with  unbounded 


strain  energy. 


SquaraJte  Arrays 

Square  fiber  arrays  coi  respond  to  the  configuration  shown  in  Figure  6, 
except  that  the  fibers  are  centered  at  points |^,  ^jand  |-^/  The 

displacement  in  regions'  A  and  a"  is  given  by  equations  (10a)  and  (10c) 

fib  b 

upon  replacing  2  by  2  . 


In  addition  vanishes  region  B.  The  singular  nature  of 


aw 


m 


ay 


and 


the  unboundedness  of  the  strain-energy  follow  in  the  same  maimer  as  in 
hexagonal  arrays  with  Vf  <  0.25. 

Conclusions 


It  has  been  shown  that  the  replacement  of  hexagonal  or  square  arrayed 
fiber  geometries  by  an  equivalent,  plate-like,  layered  geometry  results  in  a 
displacement  field  which  is  inappropriate  to  determine  the  compressive 
strength  of  fiber-reinforced  composites.  A  model  based  upon  a  plate-like 
geometry  obscures  the  three  dimensional  features  of  the  displacement  field 
and  may  give  highly  unrealistic  predictions. 
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Fig.  1.  A  fiber  reinforced  composite  medelled  as  a  two  dimensional  lamellar  region 
consisting  of  fiber  and  matrix  plates. 


Fig.  2.  Composite  compressed  in  fiber  direction..  Deflection  shape  of  all  fibers  are 
assumed  identical.  Displacement  of  matrix,  w^  varies  linearly  with  respect  to  x. 
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Fig.  3.  Hexagonally  arrayed  fibers  in  composite.  Part  in  dotted  rectangle  is  repre¬ 
sentative  volume  element. 


y 


Fig.  4.  Represcnutivc  volume  element  for  hexagonally  arrayed  composite  with  > 
0.25. 
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(b) 

Fig.  5.  Out-of*plane  displacement,  w,  of  (a)  RVE  in  Fig.  4  under  compression  (Pi 
indicates  location  corresponding  to  Fig.  4),  and  (b)  two  dimensional 
model  in  Fig.  2. 
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